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We investigate the effect of matter density fluctuations on supernova collective neutrino flavor 
oscillations. In particular, we use full multi-angle, 3-flavor, self-consistent simulations of the evo- 
lution of the neutrino flavor field in the envelope of an O-Ne-Mg core collapse supernova at shock 
break-out (neutrino neutronization burst) to study the effect of the matter density "bump" left by 
the He-burning shell. We find a seemingly counterintuitive increase in the overall v c survival prob- 
ability created by this matter density feature. We discuss this behavior in terms of the interplay 
between the matter density profile and neutrino collective effects. While our results give new in- 
sights into this interplay, they also suggest an immediate consequence for supernova neutrino burst 
detection: it will be difficult to use a burst signal to extract information on fossil burning shells or 
other fluctuations of this scale in the matter density profile. Consistent with previous studies, our 
results also show that the interplay of neutrino self-coupling and matter fluctuation could cause a 
significant increase in the v E survival probability at very low energy. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Core collapse supernovae are fantastic engines for the 
creation of large neutrino fluxes. In turn, these fluxes can 
engender large scale, collective neutrino flavor oscillations 
deep in the supernova envelope (see Refs. [THE] and see 
the review in Ref. [E] and references therein). In this 
paper we investigate a puzzling aspect of collective neu- 
trino flavor transformation in supernovae: in some cases 
a matter density fluctuation can increase neutrino flavor 
transformation rather than decrease it as simplistic mod- 
els including neutrino self-coupling seemingly predict. 

Recent numerical simulations of neutrino flavor trans- 
formation in O-Ne-Mg core collapse supernovae [17] illus- 
trate this conundrum. The results of these calculations 
agreed phenomenologically with the standard theoretical 
frame work. However, the final neutrino flavor distri- 
bution in these simulations revealed that in the normal 
neutrino mass hierarchy, neutrinos initially in mass state 
3 (the heaviest mass eigenstate) hopped to lower mass 
eigenstates with much lower probability than would be 
predicted by simple theoretical models [TU [IS] . 

In this paper we analyze the results of multi-angle, self- 
consistent, 3-flavor simulations of neutrino flavor evolu- 
tion in the neutronization neutrino burst of an O-Ne-Mg 
supernova. Here we study the particular case of flavor 
evolution of a pulse of primarily electron flavor neutrinos 
with an average energy of 11 MeV and a peak luminosity 
of 1 x 10 53 ergs _1 for a specified set of neutrino mixing 
parameters and emission spectra. 

To study why theoretical predictions of the neutrino 
mass state hopping rate differ from what is observed in 
our simulations we chose to vary only a single parameter 
in our model of the O-Ne-Mg supernova, the matter den- 



sity profile. This affords us an opportunity to conduct 
an interesting side investigation. We explore the possi- 
bility that the neutrino signal from this model could be 
used to detect features in the matter density profile of 
the supernova, assuming a knowledge of neutrino mixing 
parameters. 

Terrestrial experiments, like the proposed long baseline 
neutrino experiments, hold great promise for revealing 
key neutrino flavor mixing parameters, such as the value 
of 6*i3 and the neutrino mass hierarchy. These as yet 
unmeasured quantities influence how neutrinos change 
their flavors in the core collapse supernova environment. 

If experiments can reveal neutrino mixing parameters, 
it stands to reason that the signal from a supernova could 
be used as a probe of supernova physics. There is a rich 
physical interplay between the hydrodynamic motion and 
nuclear abundances in a supernova and the neutrino flux 
streaming out from the proto-neutron star at it's heart. 
Armed with a refined understanding of neutrino flavor 
transformation physics, it is reasonable to ask whether 
the supernova neutrino signal could be used as a probe 
of the matter density profile in a supernova at times and 
depths that are impossible to measure with optical ob- 
servations. 

For the purposes of this study we have chosen the 
following neutrino mixing parameters: neutrino mass 
squared differences Atoq = 7.6 x 10~ 5 eV 2 , Am 2 tm = 
2.4 x 10 _3 eV 2 ; vacuum mixing angles 6\2 = 0.59, #23 = 
7r/4, 6*i3 = 0.1; and CP-violating phase (5 = 0. Here we 
will concentrate on the normal neutrino mass hierarchy. 
Along with this we model the neutronization neutrino 
pulse to be a pure electron neutrino Fermi-Dirac spec- 
trum with average energy (E v ) = HMcV, a degeneracy 
parameter 7/ = 3, and luminosity L = 1.0 x 10 53 ergs -1 . 
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Previous simulations [T7] and semi- analytic work [H] 
[T5] agree broadly on the theoretical framework that 
should describe the flavor evolution of neutrinos in this 
case. In section II we discuss neutrino flavor transforma- 
tion with and without a matter density fluctuation. In 
section III we discuss the methodology of our numerical 
calculations, while in section IV we outline a theoreti- 
cal framework for collective neutrino oscillations in this 
regime. We give an analysis of our numerical results in 
section V and conclusions in section VI. 



II. NEUTRINO FLAVOR TRANSFORMATION 
WITH AND WITHOUT MATTER 
FLUCTUATIONS: A CASE STUDY 

For our particular model, neutrinos emerging from 
the neutrinosphere initially are in pure electron flavor 
states. As these neutrinos stream outward through the 
envelope, a collective effect known as the "Neutrino En- 
hanced MSW" effect (not to be confused with the MSW, 
Mikheyev-Smirnov-Wolfenstein effect) can produce mass 
state hopping of neutrinos out of the heaviest mass eigen- 
state (for both neutrino mass hierarchies) when matter 
densities are large. Nominally, the hopping rate is set 
by a comparison of the scale height of the matter den- 
sity and a characteristic neutrino oscillation length in 
the resonance region. When matter densities fall further, 
a second collective effect called the "Regular Precession 
Mode" begins. All neutrinos in this mode begin to rotate 
around an effective field in flavor space at the same fre- 
quency, regardless of their energy. Because of the v e only 
emission of the neutronization neutrino pulse, this pro- 
cess conserves the total number of neutrinos occupying 
each mass eigenstate and produces the distinctive "Fla- 
vor Swaps" or "Spectral Swaps" seen in the final neutrino 
spectra. 

This last point is extremely important. By conserving 
the number of neutrinos in each mass state, the flavor 
swaps freeze the flavor evolution history of the neutrinos 
into the final spectrum with a signature that stands out 
dramatically to an observer here on Earth. This suggests 
that such an observer might be able to simply measure 
the swap energies in a detected supernova neutrino signal 
and work backward to construct an in-situ measurement 
of the matter density profile at high densities. 

Because the neutronization neutrino pulse of an O- 
Ne-Mg supernova has been well studied and is a rela- 
tively simple case of flavor swap formation, it serves as a 
good test case to study our ability to extract information 
about the supernova envelope from a detected neutrino 
burst signal. For the inverted neutrino mass hierarchy, 
the sequence of events of neutrino flavor transformation 
produce only a single swap because only one mass state 
level crossing is present [HI [TSJ [17] . For the normal neu- 
trino mass hierarchy, two swaps are produced because 
two separate level crossing populate all three mass eigen- 
states with neutrinos (for small $13 this can be reduced 



to a single swap via the complete depopulation of mass 
state 3) [HHHl HI]. 

The density profile used for this study comes from a set 
of simulations by Nomoto [TH1 US] ■ This profile is typical 
of what a mid-collapse O-Ne-Mg supernova might pro- 
duce at the epoch of the neutronization neutrino pulse, 
r~> 10 ms post bounce. A feature that this profile pos- 
sesses is a small bump in the matter density which is 
created by the star's helium burning shell in a range of 
radius bounded by r ^ 1080 km — 1100 km . This fea- 
ture is known to cause neutrinos at low energies to pass 
through multiple MSW resonances at the Am^tm mass 
scale, and has been discussed in [T3 | H5 | H7]. 

To test our ability to detect a simple feature such as 
this, we conducted a pair of simulations. The first with 
the original density profile, called "Bump" , and the sec- 
ond with a synthetic density profile where the He burning 
shell feature has been removed, called "No Bump" . Fig- 
ure [T] shows the electron number densities with these two 
profiles plotted side by side. For the neutrino-electron 
forward scattering potential H c (hereafter the "matter" 
potential) the associated scale height at resonance is, 



Figure [2] shows the matter potential scale height for both 
profiles, evaluated at the MSW resonance location for 
each neutrino energy bin. 

Figures [3] and [4] show the results of our calculations for 
the flavor transformation of electron neutrinos emitted 
during the neutronization neutrino burst. Figure [3] shows 
the results of the original, Bump, density profile. Figure[4] 
shows the results of a simulation using the No Bump 
density profile. 

The aim of the second simulation, with the No Bump 
density profile, was to study whether we could detect a 
signature from features of the matter density profile us- 
ing the neutrino flavor transformation signal. Post pro- 
cessing of this data led to a surprise. The total num- 
ber of neutrinos that remain in the heavy mass eigen- 
state decreases when the bump in the density profile is 
removed. Explicitly, the heavy mass eigenstate (mass 
state 3) survival probability, Pjj, for the two cases are 
P H Bump = 0.852 and Bump = 0.759. This is a coun- 
terintuitive result. We expected that the removal of the 
bump from the original density profile would have pro- 
duced flavor evolution that was more adiabatic, leading 
to a greater survival probability for the No Bump profile. 

The effect that the removal of the helium burning shell 
has on the v c survival probability is shown in Figure [5] 
Note in this figure that there is an enhanced survival 
probability for 14's at low energy for the Bump profile, 
and that the flavor swap energy is lowered for neutrinos 
in the No Bump profile. 

Collective flavor evolution prior to the onset of the reg- 
ular precession mode was thought to be described by the 
Synchronous MSW effect, as opposed to the Neutrino 
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FIG. 1: The electron number density for two matter profiles 
plotted as functions of radius in the resonance region for the 
Arria tm mass state splitting. The solid line indicates the orig- 
inal matter density profile of Refs. [181119] , called Bump. The 
dashed-dotted line indicates the artificial density profile with 
the bump artificially removed, called No Bump. 
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FIG. 2: The scale height of the neutrino-electron forward scat- 
tering potential evaluated at the MSW resonance location for 
each neutrino energy. The solid line indicates the first MSW 
resonance scale height of neutrinos moving through the orig- 
inal matter density profile of Refs. [181 119] , while the dashed 
line indicates the scale heights of the multiple resonances. The 
dash-dotted line indicates the MSW resonance scale heights 
of neutrinos moving through the No Bump density profile. 



Enhanced MSW effect [H Q]3 Q7J . While there is a tech- 
nical difference between these two modes of neutrino fla- 
vor transformation, the neutrino flavor transformation 
survival probabilities (and consequently the Swap ener- 
gies) are quite similar for both cases, which led to the 
initial confusion. In both cases the resultant value of 



Ph is determined by the evolution of a single represen- 
tative neutrino flavor state, which is the flavor state of 
the collective ensemble of neutrinos. Ideally, neutrinos 
remain so closely aligned with this collective flavor state 
that they do not "feel" the neutrino self coupling poten- 
tial and, as a result, they behave as a single neutrino 
experiencing the matter-driven MSW effect. 



III. METHODOLOGY 

The simulation results analyzed here were produced 
with two numerical codes used for simulating neutrino 
flavor evolution. These codes, the FLAT code and the 
BULB code, and related schemes to solve for the flavor 
evolution of core collapse supernova neutrinos, are dis- 
cussed in Ref. [20] , 

In order to parallelize the nonlinearly-coupled differ- 
ential equations which describe neutrino flavor evolution, 
BULB employs a spherically symmetric representation of 
the region above the neutrino sphere. All neutrinos are 
assumed to be emitted from a hard spherical shell, and 
propagate through a one dimensional distribution of mat- 
ter. These choices allow the neutrino emission to be bro- 
ken down and grouped by species, energy, and emission 
angle at the neutrino sphere. Here we define the emission 
angle to be the angle between the neutrino direction and 
the vector normal to the surface of the neutrino sphere at 
the neutrino emission point. To initialize the simulation, 
neutrinos are allocated to each energy-angle bin accord- 
ing to the species-specific luminosity and neutrino energy 
spectrum characteristics. From there, BULB employs a 
second order predictive-corrective algorithm to compute 
the flavor evolution of the neutrino states. In order to 
check for convergence, a second round of computations 
are made with a step size At/2. The final neutrino flavor 
states from the double iteration are then compared to the 
final states of the original step to verify that they agree 
to within a predefined error tolerance, usually chosen to 
be 1 part in 10 s . Convergence of the overall calculation 
is checked by comparing results with different error tol- 
erances and differing numbers of energy and angle bins. 

We have validated our simulations by performing them 
with both codes and comparing our results [171120]. Using 
the same set of initial conditions, both codes agree with 
each other at the level of 0.1% when comparing the final 
neutrino flavor states. For a more detailed description 
of the inner workings of the BULB and FLAT codes, see 

Refs. quehch]- 



IV. THEORY 

For our simulations we assume a pure v e burst emitted 
from the neutrino sphere ar R v — 60 km with a total lu 
minosity of L v = 10 53 



erg s 1 and a normalized spectrum 
1 E 2 



F 2 fa,,) 1* exp {E/T v - Vu ) + 1 ' 



(2) 
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FIG. 3: Bump. Left panel: electron neutrino survival probability P Vb v b (color/shading key at top left) for the normal mass 
hierarchy is shown as a function of cosine of emission angle, cosi?r, and neutrino energy, E in MeV, plotted at a radius of 
r — 5000 km. Right: mass basis (key top right, inset) emission angle-averaged neutrino energy distribution functions versus 
neutrino energy, E. The dashed curve gives the initial u e emission angle- averaged energy spectrum. A kink in the density 
profile used, taken from Refs. |T5] 119] , leads to multiple MSW resonances for low energy neutrinos. 




FIG. 4: No Bump. Left panel: electron neutrino survival probability P Vt , Vc (color/shading key at top left) for the normal mass 
hierarchy is shown as a function of cosine of emission angle, cosi?r,, and neutrino energy, E in MeV, plotted at a radius of 
r — 5000 km. Right: mass basis (key top right, inset) emission angle-averaged neutrino energy distribution functions versus 
neutrino energy, E. The dashed curve gives the initial u c emission angle-averaged energy spectrum. The kink in the density 
profile taken from Refs. [181 119j . has been artificially removed from the density profile used in this simulation. 



where r\ v = 3 and T u = 2.75 MeV. This corresponds to an 
average v energy (E u ) = F 3 (??„) T„/F 2 (?7„) = 11 MeV. 
Here 



exp (a; —rj u ) + l 



dx 



(3) 



In the single-angle approximation, the effective total neu- 
trino number density at r > R v is 



n v (r) = 



i-vi-wr 



L V R V 
167r(^)r 4 ' 
(4) 
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FIG. 5: The electron neutrino survival probability as a func- 
tion of neutrino energy shown at the final radius r = 5000 km. 
The solid line indicates the survival probability for neutrinos 
after moving through the Bump profile, while the dashed line 
indicates the survival probability for neutrinos after moving 
through the No Bump profile. 



where the approximate equality holds for r 3> R v , and 
where we set h = c = 1. 

Because the salient features of our results are confined 
to the Sm 2 tm mass squared mixing scale, we will confine 
the following discussion to a two neutrino flavor mixing 
scheme. Following the convention of [22], we take each 
neutrino flavor state with energy E v and represent it as 
a three dimensional neutrino flavor isospin (NFIS), 



e£/2, for u e , 
-el/2, for v x , 



(5) 



where e l z is the unit vector in the z-direction for the neu- 
trino flavor basis and oj is the vacuum oscillation fre- 
quency u — 8m 2 atm /2E y . We focus on the normal mass 
hierarchy and we take the effective vacuum mixing angle 
6 V rj O13 = 0.1. The evolution of a NFIS s w is governed 
by 



dr 



wH v + H e - M (r) / s,,/ dE* 
Jo 



(6) 



where H v = cos20 v e£ — sm29 v e f x , H G = — ^/2Gpn e (r) e*, 
fi(r) — 2^2G F n v (r), and E u > = £m 2 tm /2a/. For conve- 
nience, we define 



g{u) = S J^ fviEu) 



(7) 



and 



It follows from Eqns. that 
d 



dr 



-S = 



ujg (oj) s^dui x H v + S x H c 



(9) 



As g (oj) is concentrated in a finite range of oj, to ze- 
roth order we approximate g (oj) ~ 5 (oj — (oj)), where 
(lu) = J ujg (uj) duj. This is a fair approximation for the 
particular case we treat, namely, low energy ^ e 's in the 
neutronization burst ((E u ) — 11 MeV). With this ap- 
proximation, the zeroth-order mean field S(°) is defined 
through 



dr 



S (0) = g(0) x faajj + H j = g(0) x jj 



MSW 



(10) 



The evolution of S*- -* is the behavior of the system in 
the high luminosity, "Synchronized", limit. (In the Syn- 
chronous MSW effect, all s w are aligned with S, and orbit 
around it. Note, this idealized situation does not occur in 
the presence of matter |23j.) For the neutronization neu- 
trino burst luminosity and matter density profile we use, 
our calculations take place below this luminosity regime. 
We do not observe individual s w orbiting S [24], but we 
do observe that individual s w remain closely aligned to S 
andSW. 

The evolution of S (0) is the same as that of an idealized 
v e with E v — 6rn^ tm /2(oj) = 8.53 MeV undergoing the 
usual MSW effect. This is the collective NFIS that all the 
the neutrinos will follow during the Neutrino Enhanced 
MSW effect. Now we can approximately solve for the 
evolution of Sm from 



d 

— < 

dr 



wH v + H c 



/' 



(r) S ( °) 



(11) 



We can use the solution of Equation [10] to obtain the 
first-order mean field from the definition of S in the 
first expression in Equation [HJ and then re- calculate the 
evolution of s w from the above Equation 10 but with S' 1 ) 
replacing S^ -*. This procedure can be repeated until the 
results converge. 

Of course, the procedure outlined above is not recom- 
mended as a numerical method, but instead points to an 
analytic approach to understand collective oscillations. 
Based on the MSW effect, S*- ^ goes through resonance 
at n c w 1.09 x 10 27 cm -3 corresponding to r « 1100 km in 
both simulations. For simplicity, let us consider a smooth 
distribution for n Q (r). Before the resonance, the evolu- 
tion of S*- -* is somewhat adiabatic and we can take 
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■- (cos20 n 



sin 29„,e 



(12) 
(13) 



where 



cos 29„ 



cos 26 v - y/2G F n e 



S= / s u f v (E u )dE u 



s^g (oj) dcu . (8) \J ((oj) cos26» v - V2G F n c ) 2 + ((oj) sin20 v ) 2 



(14) 



6 



10" 



03 i n- 
CO 



10" 



-17 

1000 



Radial 



v^GWc + ^r) (S-e*) 
£ ros = 3.2 Me V 
E ies = 4.2 MeV 



1020 1040 



1060 



1080 



Radius (km) 




1100 



10" 



10" 



10" 



10" 



10" 



10" 



1000 



Tangential 



V2G F n c + ii (r) (S ■ e { z ) 
E ICS = 2.0 MeV 
E Tes = 2.5 MeV 




1020 



1040 



1060 



1080 



Radius (km) 



1100 



FIG. 6: Left panel: the combined matter and neutrino self-coupling forward scattering potentials experienced by radially 
emitted neutrinos as a function of radius. Right: the combined matter and neutrino self-coupling forward scattering poten- 
tials experienced by tangentially emitted neutrinos as a function of radius. The dotted and dashed-dotted lines indicate the 
upper and lower neutrino energies, respectively, to experience multiple neutrino-background enhanced MSW resonances in the 
desynchronized limit. 



sin 26 m = The evolution of s w at densities higher than the resonance 

(lj) sin 28 density for S 1 - ' is then governed by 

(15) 



^ ((lj) cos2(9 v - V2G F n c ) 2 + ((lj) sin20 v ) 2 



dr 



(lj cos 26 v — V2Gpn c 



~ cos 26> 



.) 4 - ( wsi 



sin 29 v + — sin 29 n 



x 



(16) 
(17) 



This equation also defines H w . Note that the NFIS where 
evolution described by this equation is similar to the 
usual MSW effect but with modified diagonal and off- 
diagonal terms. Note especially that the off-diagonal 
term (jti/2) sin 28 m = /i/2 is large at the resonance lo- 
cation for S^ ). If s u evolves adiabatically, then 

" ' (18) 



2|H 
1 



( v cos20 w e* -sin 26^4) , ( 19 ) 
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cos 26 u 



(uj) cos 29 v - y/2G F n c + (/i/2) cos 29„ 



\J ((uj) cos26> v - y/2G F n c + (/i/2) cos20 m ) 2 + ((us) sin26» v + (/i/2) sin26> m ) 2 



(20) 



sin 26* 



(ui) sin26> v + (/i/2) sin26! n 



^ ((w) cos2# v - v^Gf^c + (/i/2) cos26» m ) 2 + ((w) sin20 v + (/i/2) sin26> m ) 2 



r 



(21) 



If s w goes through resonance non-adiabatically at den- 
sities above the resonance density for S^°\ then subse- 
quently s u will no longer stay anti-aligned with H w . This 
change in alignment corresponds to neutrinos jumping 
between mass states. Confirmation of this simple picture 
is borne out by our numerical simulations as can be seen 
in Figures [3] and |4j where the Bump profile exhibits a 
population of low energy mass state 2 (i/ 2 ) neutrinos. In 
the Bump density profile, Figure |3j the helium burning 
shell produces multiple MSW-like resonances for low en- 
ergy neutrinos. These multiple resonances, illustrated in 
Figure |6j are non-adiabatic, engendering further loss of 
alignment. This leads to a population of low energy neu- 
trinos occupying mass state 2, which do not recover their 
alignment with H w . 

By way of contrast, in Figure [4] there is no population 
of low energy neutrinos in mass state 2 because the ab- 
sence of the helium burning shell in the No Bump profile 
makes the evolution of these neutrinos adiabatic, hence 
they remain aligned with H w . In both cases, NIFS's 
over a wide range of w experience significant evolution at 
densities higher than the resonance density of S(°K 

If S'°) goes through resonance adiabaticaUy, the above 
description of NFIS evolution can be extended to lower 
densities. Note that cos2# m changes from « — 1 at high 
density to at resonance and to sa cos 29 v at low density, 
consistent with a simple MSW picture. The diagonal 



termv2GF^ c — (/i/2) cos 2(9 m becomes at some radius 
and all s u go through resonance before this radius. 

However, if goes through resonance non- 

adiabatically, the situation becomes more complicated. 
For illustration consider the regime of low n c , lower than 
the S^ ) resonance density. The non-adiabatic evolution 
of S^ * 1 means that it no longer stays anti-aligned with 
Hmsw ~ (w)H v in this regime. Instead, we have 



S (0) H v = icosa= 1} 



hop 



(22) 



where a is the angle between S^ ^ and the external field 
H v , and Phop is the probability for S' ) to hop from be- 
ing anti-aligned before resonance to being aligned with 
Hmsw after resonance. We can take 



S (0) 



Pi 



hop 



H v 



(23) 



where S? is the component perpendicular to H v with a 

(o) 



magnitude « yj (1/2) - [P hop - (1/2)] . Note that S 
precesses around H v with an angular velocity — (u;)H v . 

In a frame co-precessing with S^°\ the evolution of 
is governed by 



dr 



Suj X { 

s.„ x H! :< " 



- (w) -fi(r) ^P ho p - 



H v - M (r)sf } 



(24) 

(25) 
(26) 



The evolution of s w is expected to be adiabatic and, as 
a result, the angle between s w and H^°~ prc stays fixed. 
This angle depends on the relative directions of s u and 
H co- P rc right after the resonance f The latter di- 
rection depends on the exact direction of S_l coming out 
of the resonance. At large radii, H™ _prc simply becomes 
(ui — (uj) ) H v . Those with uj < (uj) that are approxi- 



mately aligned with prc right after the resonance of 
S 1 - -* are nearly fully converted into u x . 

Neutrinos participating in the Neutrino Enhanced 
MSW effect have NIFS's which are closely aligned with 
the collective field S, so their flavor evolution will be a 
close match to that of [23 . Specifically, they will 
co-precess with the effective field, even when the flavor 



evolution of and S are non-adiabatic. Examples of 
this co-precession can be viewed in movies which are on 
our website [23] . 

In this co-precession picture the probability to re- 
main in the heavy mass eigenstate, Ph, should depend 
on the three quantities relevant to the MSW evolu- 
tion of S' ); vacuum mixing angle 9 V ; collective os- 
cillation frequency (w); and the matter potential scale 
height i?H at the location where is at resonance, 
(u>h) cos 28 v = 2v / 2G < Fn c (r). The first and second quan- 
tities are identical for both of these simulations. Only 
i?H changes when the burning shell feature is removed, 
with P H ump = 25.3 km and Bump = 28.5 km. 

Following the evolution of S' ' through the envelope 
of matter around the proto-neutron star, the probabil- 
ity of a neutrino with oj = (uS) to hop out of the heavy 
mass eigenstate will be given by the double exponen- 
tial Landau-Zener hopping probability, Ph op = (1 — Ph)> 
with: 



Pi 



hop 



exp (27rP H ( w ) cos2 #v) — 1 
exp(27rP H (w)) - 1 ' 



(27) 



Given that the critical scale height of the matter profile is 
slightly smaller for the Bump profile, the above equation 
implies that there should be a larger probability to hop 
out of mass state 3 in the presence of the helium burning 
shell. This is why the total number of neutrinos remain- 
ing in mass state 3 naively is expected to increase for the 
No Bump simulation. Equation |27| yields a prediction 
that for the Bump density profile P£ rcdictod = 0.68, and 
for the No Bump density profile predicted = 0.71. 

Using the final emission spectra from Figures [3] and [4] 
to calculate Ph for both simulations, as we would in an- 
alyzing an actual supernova neutrino signal, produces 
somewhat different results. For the No Bump calcula- 
tion, pObscrvcd _ q 7g Thig in di ca tes a slightly more 
adiabatic than we had calculated above, and would lead 
an observer to deduce a larger matter scale height in 
the collapsing core of the supernova than was actu- 
ally present. More strikingly, the Bump profile exhibits 
pObscrvcd _ q gg^ implying fully 17% more neutrinos re- 
main in mass state 3 than predicted. This would lead 
an observer interested in the envelope to grossly miss- 
calculate the electron density scale height, arriving at a 
number nearly twice the actual value. We will endeavor 
to understand why it is that our model of neutrino fla- 
vor transformation seems to have led us astray when at- 
tempting to work backward from our observed signal to 
the matter density profile of the collapsing star. 



V. ANALYSIS 

Though the simple neutrino transformation model pre- 
sented above is successful in many respects, we do not at- 
tempt to provide an exhaustive proof of the model, only 
to point out that it offers a straightforward explanation 



of the puzzling aspects of our result. Furthermore, our 
numerical calculations show that one must be careful in 
applying the theory of the Neutrino Enhanced MSW ef- 
fect to infer information about the envelope of the proto- 
neutron star. 

First, even in early (e.g. Ref. [7J) treatments of neu- 
trino flavor evolution it was evident that flavor diago- 
nal neutrino-neutrino forward scattering potential, B (r), 
would alter the position of the MSW resonance position 
for a given neutrino energy E v and mass splitting Am 2 
because the resonance condition is 



Am 
2E 



cos26» v = A(r) + B{r) , 



(28) 



where A (r) = ^/2Gfti c (r) is the matter potential at ra- 
dius r. Neutrino propagation through the MSW reso- 
nance at the shifted position in general will result in an 
altered survival probability because the scale height of 
the combined potential will be different. However, this 
simplistic analysis is completely inadequate because it is 
the neutrino-neutrino flavor off-diagonal potential which 
in part determines adiabaticity [71 125). This potential, 
in turn, is sensitive to the relative x — y plane phase 
angles of the individual neutrino NFIS's, necessitating a 
self-consistent collective oscillation treatment. 

Shown in Figures [7] and [8] is the evolution of the Neu- 
trino Enhanced MSW effect collective NFIS S and evo- 
lution of the zeroth order approximation S'°). Each 
of these are given for the Bump and No Bump pro- 
files. These figures show the opening angle between 
the collective NFIS and Hmsw- To lowest order, the 
collective NFIS S from Eqn. [8] follows the alignment 
of S*- ', starting anti-aligned with Hmsw and undergo- 
ing a mild change of alignment as the system passes 
through resonance, with the hopping probability given 
by Phop = 1/2(1 + cos a). However, there is a small 
difference between the final alignment of S' - 1 and S for 
all of the calculations. The full numerical calculations 
reveal that mass state hopping is more adiabatic than 
our zeroth order approximation would lead us to believe, 
differing by « 5 — 10 % from the hopping probability as- 
sociated with the evolution of S^°\ 

The reason that the full calculations exhibit less mass 
state hopping can be found in the imperfect alignment 
of individual s w with S. The Neutrino Enhanced MSW 
model predicts that individual s w tend to stay aligned 
with S as the collective mode passes through resonance 
and subsequently begin to orbit around H v as the system 
transitions to the regular precession mode. In Figure [9] 
we show the average opening angle between individual 

s w and S, taken to be cos6> = • Sj, for the Bump 

and No Bump profiles. On average individual s w and 
S remain aligned to within a few percent throughout the 
resonance region in both the Bump and No Bump simula- 
tions, which shows that this is indeed the correct physical 
picture. However, as we have mentioned previously, the 
individual s w pass through resonance at slightly higher 
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FIG. 7: Bump: The opening angle a between the collective 
NFIS and Hmsw for the Bump profile, plotted as a function of 
|Hmsw| / |H V | as the system moves through resonance. The 
idealized NFIS (dotted-dashed line) shows the evolution of 
S' - 1 in the absence of any neutrino self-coupling. The solid 
line and dashed line show the evolution of S as calculated in 
multi-angle and single-angle simulations respectively, includ- 
ing the neutrino self-coupling potentials. 



FIG. 9: The average alignment of individual neutrino po- 
larization vectors, s u , with the collective polarization vec- 
tor, S. The solid black line shows the "Bump" profile from 
Refs. [181 119] and the dotted red line shows the modified "No 
Bump" profile. The shaded region indicates the physical po- 
sition of the Helium burning shell density feature present in 
the Bump profile. 
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FIG. 8: No Bump: The opening angle a between the collec- 
tive NFIS and Hmsw for the No Bump profile, plotted as a 
function of |Hmsw| / |H V | as the system moves through res- 
onance. The idealized NFIS (dotted-dashed line) shows the 
evolution of S*- ' in the absence of any neutrino self-coupling. 
The solid line and dashed line show the evolution of S' ' as 
calculated in multi-angle and single-angle simulations respec- 
tively, including the neutrino self-coupling potentials. 



densities than S^ ^. The individual s u are slightly mis- 
aligned and this means that when the individual neu- 
trino states are at resonance, i.e. cos26> w = 0, the col- 
lective mean field S is not yet at resonance itself, i.e. 



From Equation 20 we see that for an individual neu- 



trino state at resonance the adiabaticity of the mass state 
hopping will not be determined entirely by the matter 
potential if cos 26* m ^ 0, which is precisely the result we 
recover from Equation 28 when A > B with B ^ 0. (By 
contrast, the flavor evolution through resonance for S' -* 
is determined entirely by the matter potential.) Indi- 
vidual s w experience some fraction of the neutrino self- 
coupling potential fi. This comes from the fact that our 
approximation g (u>) ~ S (uj — (lo)) is a rather gross ap- 
proximation. In reality, the function g (oj) has a finite 
width. However, we find that this approximation pro- 
duces results that match well with our calculations and 
individual s u track the evolution of S'* -' more closely than 
one might expect given the width of our initial v c distri- 
bution. 

A serious source of potential error in interpreting the 
swap signal for the neutralization neutrino burst comes 
from rapid fluctuations in the matter potential. We have 
shown the evolution of the system follows S'°) closely. 
In turn, S^ ' experiences only the ordinary MSW effect 
in it's flavor evolution. This is illustrated clearly by the 
similarities in the observed mass state 3 survival prob- 
abilities and the trajectories of S and S' ' in Figures [7] 
and [sj For the No Bump profile pObsorvcd = J6; and 

0.68. For the Bump profile pObscrvcd = .85, 

However, as we have mentioned in 
the previous section, this result for the Bump profile is 
in stark disagreement with the prediction of the double 
exponential Landau-Zener hopping probability of Equa- 



ls 

and P]f <0> = 0.80. 
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FIG. 10: No Bump: The precession of the collective polariza- 
tion vector, S, and the zeroth order effective field S' ', about 
the e| axis as neutrinos move through resonance in the No 
Bump profile. The solid line shows the azimuthal angle, (j>s, 
for the collective polarization vector S. The dashed-dotted 
line shows the azimuthal angle, s <o) , for the zeroth order ef- 
fective field S ° '. The dotted line shows the value of # m for 
S and S'- ' 1 for reference, with a diamond symbol located at 
9 m — 7r/4 where the system is at resonance. 



tionEH 

As has been know for some time, turbulent matter den- 
sity fluctuations can produce flavor "depolarization" for 
MSW neutrino flavor transformation [5BH33]. Broadly, 
these turbulent fluctuations can produce mass state hop- 
ping that does not agree with what one would deduce 
from Eqn. [27] using the gross scale height of the matter 
potential. 



By considering the derivative of Eqn. 10 it can be seen 
that there is a restoring force that causes S^ ^ to or- 
bit Hmsw as the sys tem evolves. Another term in the 
derivative of Eqn. 10 drives precession of about the 



e z axis any time the matter potential is changing non- 
adiabatically. 

In the traditional MSW framework, this precession is 
dependent on the energy of each neutrino, as the neutrino 
energy sets the relative size of the terms in Eqn. [l"0| s 
derivative. As a result, if the matter density profile is 
turbulent, the alignments of individual neutrino NFIS's 
can be scattered throughout the flavor space, hence the 
term "depolarization" . 

The matter driven precession figures into the hopping 



probability given by Eqn. 27 because the double expo- 
nential hopping probability is derived from the Landau- 
Zener two level hopping problem using the boundary 
condition that the precession of S^ ) originating from 
this term can be taken to be zero before the system ap- 
proaches resonance (35] [36] . If there is significant mat- 
ter driven precession for S and S^ -* prior to resonance, 



FIG. 11: Bump: The precession of the collective polarization 
vector, S, and the zeroth order effective field S' ', about the 
e z axis as neutrinos move through resonance in the Bump 
profile. The solid line shows the azimuthal angle, 4>s, for the 
collective polarization vector S. The dashed-dotted line shows 
the azimuthal angle, 4> S (o) , for the zeroth order effective field 
S (0) . The dotted line shows the value of m for S and S (0) for 
reference, with a diamond symbol located at m = 7r/4 where 
the system is at resonance. The shaded region indicates the 
physical position of the Helium burning shell density feature 
present in the Bump profile. 



the Landau-Zener hopping probability. 

These issues can be explored by examining the projec- 
tion of the collective field S (or S' ^) in the e^ — e f y plane. 
Define the angle made by this projection and e^. to be 4>. 
Figure [To] shows the value of <p in the No Bump profile for 
S and S^ u ) through the resonance region. Prior to reso- 
nance, when 9 m i: tt/2, there is no appreciable precession 
of either S or S*- -*. This is not particularly surprising, as 
the observed hopping probability of S 1 - ' matches exactly 
with Eqn. |27[ with S exhibiting slightly less hopping as 
we have discussed above. 

where we show the evolution of 6 for S 



Figure 11 



Eqn. 27 will not be the appropriate analytic solution for 



and S(°) in the Bump profile, exhibits very different phe- 
nomenology. The system starts out identically to the con- 
figuration of the neutrinos in the No Bump profile, with 
no matter driven precession before the bump is reached. 
However, the sudden increase in the local matter den- 
sity at r ~ 1086 km brought on by the helium burning 
shell drives rapid precession of S and S'°) in a counter 
clock-wise direction. At the same time, the value of 9 m is 
brought back up to the pre-resonance value of 6 m ^ it/2. 
This resets the clock, so to speak, on Hmsw an< i the 
boundary conditions for S and S( ). The collective NFIS 
is clearly precessing rapidly at the boundary of the MSW 
resonance region, meaning that our use of the double ex- 
ponential hopping formula is not appropriate. 

This phenomenon is a simplistic example of the turbu- 
lence driven flavor depolarization discussed in Ref. [551 
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1501 1551 15i] , where we are following the evolution of a sin- 
gle NFIS, S(°\ as it is scattered to a random position in 
flavor space. What is surprising about our calculation is 
that we have observed that the collective Neutrino En- 
hanced MSW effect persists and maintains coherence in 
multi-angle calculations, in spite of rough handling by 
the matter potential. The collective NFIS, S (and con- 
sequently the individual neutrino s w 's), tracks the evo- 
lution of S' ) and achieves a final alignment in flavor 
space that is not relatable to the simplistic prediction of 
Eqn. [27] Our calculations suggest that the Neutrino En- 
hanced MSW effect is susceptible to the same turbulence 
driven precession as the ordinary MSW effect. 

This is unfortunate from the perspective of attempts 
to use the swap signal to probe the matter density pro- 
file. In the particular case of the neutronization neutrino 
pulse, the resultant swap energy may or may not be re- 
latable to the matter density profile, depending on the 
"smoothness" of the actual matter density profile of the 
collapsing star. In the case of the Bump profile, our 
observed survival probability of p® ump = 0.852 would 
lead an observer to infer a matter density scale height of 
^Bump, infoiicd _ 45 compared to the actual scale 

height of i?^ ump = 25.3 km, an 85% error. In principle, 
this fluctuation driven precession might alter the survival 
probability Ph and, consequently, push the swap energy 
up or down depending on the particular realization of the 
density fluctuations [2"51 |501 |551 151] . 

In the example we study here, the low energy neutrinos 
experiencing multiple MSW resonances would be almost 
invisible to neutrino detectors designed to collect super- 
nova neutrinos, since these detectors might have low en- 
ergy thresholds of ~ 5—10 McV. However these neutrinos 
provide our only clue about matter density fluctuations. 
While we have been able to use our detailed knowledge 
of the neutrino flavor states as they evolved through the 
resonance region to correctly deduce what has transpired, 
an observer on the Earth would not have access to such 
privileged information. 

VI. CONCLUSION 

We have elucidated the rich interplay of collective neu- 
trino oscillations and the underlying matter density pro- 
file in a supernova. Our calculations reveal a heretofore 
unrecognized aspect of collective oscillations, namely the 
matter density fluctuation-driven precession of the col- 
lective ensemble of neutrino NFIS's (e.g. please view 
the movies located at [24]). The implications of this for 
neutrino-detection astrophysics are somewhat negative. 

In particular, we have shown that during the neutron- 
ization neutrino pulse epoch of an O-Ne-Mg core-collapse 
supernova there are strong limitations on an observer's 
ability to use the swap signal to probe fluctuations in 



the matter density profile of the collapsing star (e.g. the 
density ledges left by fossil burning shells). If the mat- 
ter density profile of the progenitor O-Ne-Mg star is not 
smooth, reverse engineering the swap energy to find a 
matter density scale height may not necessarily be pos- 
sible. 

The problem hinges on an observer's ability to detect 
neutrinos which have energies low enough that their fla- 
vor evolution histories would be disconnected from collec- 
tive flavor oscillations in the presence of matter density 
fluctuations. If low energy features such as the one seen 
in Figure [3] are observed during the neutronization neu- 
trino pulse epoch, an observer could at the very least de- 
duce the presence of a fluctuating matter density profile. 
However, this requires either a neutronization neutrino 
pulse that has average neutrino oscillation frequency, (lo) , 
larger than expected from progenitor models, or neutrino 
detectors that have neutrino energy detection thresholds 
below 5 MeV. 

The model presented here suggests further investiga- 
tion. It is likely that the Neutrino Enhanced MSW effect 
could provide a secondary probe of the neutrino lumi- 
nosity during the neutronization burst. Furthermore, v e 
fluxes during this epoch with (u) different from the value 
that we have used here may also produce unique flavor 
transformation signatures that could potentially be used 
as an independent measure of the u e temperature. 

It should be noted that these effects are limited to su- 
pernova neutrino systems that experience the Neutrino 
Enhanced MSW effect, which is not universal to all su- 
pernova neutrino systems. As the neutrino flux from the 
proto-neutron star evolves past the neutronization neu- 
trino burst phase, entirely different collective effects are 
capable of producing swaps which may not be susceptible 
to turbulent matter density fluctuations. 

Nevertheless, the study presented here gives new in- 
sights into collective neutrino oscillations in supernovae. 
This and other neutrino flavor oscillation studies may 
help drive synergy between the laboratory neutrino 
physics enterprises and observational and theoretical as- 
trophysics. 
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